On the 2-point function of the O(N) model 

Federico BenitezQ 
Instituto de Fisica, Facultad de Ciencias, 
Universidad de la Republica, Igud 4225, 11000, Montevideo, Uruguay 

Ramon Mendez Galanii and Nicolas WscheboJi] 
Instituto de Fisica, Facultad de Ingenieria, Univ. de la Republica, 
J.H.y Reissig 565, 11000 Montevideo, Uruguay 
(Dated: February 1, 2008) 

Abstract 

The self-energy of the critical 3-dimensional O(N) model is calculated. The analysis is per- 
formed in the context of the Non-Perturbative Renormalization Group, by exploiting an approx- 
imation which takes into account contributions of an infinite number of vertices. A very simple 
calculation yields the 2-point function in the whole range of momenta, from the UV Gaussian 
regime to the scaling one. Results are in good agreement with best estimates in the literature for 
any value of N in all momenta regimes. This encourages the use of this simple approximation 
procedure to calculate correlation functions at finite momenta in other physical situations. 
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I. INTRODUCTION 



The O(N) scalar model describes many phenomena in a wide range of physical sit- 
uations. Besides the N — 1 case, which corresponds to Ising like systems, with a wide 
range of applications as, e.g., liquid-gas transition, the N = 2 model describes superfluid 
Helium, N = 3 can be used to study ferromagnets, N = 4 allows the study of the Higgs 
sector of the Standard Model at finite temperature, and the N = case describes the 
physics of some polymers As a natural consequence, a huge amount of work has been 
devoted to the study of this family of models. Leaving aside the d = 2 case, where specific 
methods exist, most of the existing results correspond to thermodynamical properties, as 
critical exponents or phase diagrams, i.e., physical quantities encoded in correlators at 
small external momenta, e.g., the effective potential. With this goal, very complicated 
techniques, as resummedperturbative calculations carried up to 7-th order [l|, Q , high- 
tetnperatnte expansions Bflfl, o r Monte-Co .netfiods QflflQg , were used. When 
instead trying to get ph ysl ea, quantities depending on finite —a, snch as the self- 
energy of the model, fewer results can be found in the literature (see for example [9[ and 
references therein) . All these calculations suffer from a common difficulty, which is general 
to a vast class of problems: it is extremely nontrivial to deal with systems having highly 
correlated components. In this sense, as the O(N) model is simpler than most other such 
problems, it has been largely used as a testing ground for the development of calculation 
schemes in non-perturbative contexts. 



The Non-Perturbative Renormalization Group (NPRG) [10, QjJ, ll2j, [13|, [14j is a general 
framework conceived to deal with this kind of situations. It is based in an infinite set of 
exact equations giving all renormalized correlation functions between the various compo- 
nents of a given system. Naturally, as one has to deal with an infinite tower of coupled 
differential equations, in order to solve them the use of approximations is unavoidable. 
Several years ago, a systematic approximation scheme was developed Q, which 
allows for the solution of this set of equations in a particular case: The so called derivative 
expansion (DE) is based in an expansion in the powers of the derivatives of the fields. 
Even if there is no formal proof of its convergence, the DE has provided very competi- 
tive results in problems where only small (eventually zero) external momenta play a role. 
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Among many other applications (see, e.g. 
the 0(N) model 16, 3, 3, 21 



scheme in the N = 1 Ising case 



22 



13, 



181]). the approximation was applied to 



231 ] . even up to the next-next-to-leading order of the 



24[ , yielding at this order critical exponents of a similar 
quality as those obtained using 7-loops resummed perturbative calculations. 

On the other hand, when trying to describe phenomena involving all modes, the 
situation is different. For example, to get the transition temperature of a dilute gas to 
a Bose-Einstein condensate, one needs the self-energy of the 0(2) model in 3-dimensions 
at arbitrary momenta 25| ; relevant information comes from the intermediate momentum 
region between the IR and the UV ranges. Within the NPRG, only calculations includin. 
a finite number of vertices 26| had been considered up to now, either in O(N) 27|, [28 



29, 



30. 



3jJ or in more involved problems such as QCD 32l. l33|. 



Recently, a general approximation scheme suitable to get any n-point function at any 



finite momenta within the NPRG has been proposed [34J. The strategy has many inter- 
esting similarities with DE. First, it can be applied in principle to any model. Second, 
although the approximation is not controlled by a small parameter, it can be systemati- 
cally improved. Furthermore, this strategy reproduces both perturbative and DE results, 
in their corresponding limits; for example, if solving the 2-point function flow equation at 
the leading order (LO) of the procedure, the 2-point function includes all 1-loop contribu- 
tions while the effective potential includes all 2-loops ones. It is possible to apply, on top 



of the approximation presented in 



34[ , an expansion in powers of the field, as freque ntly 



done in DE: at least in the studied case, this expansion seems to converge rapidly 35]. 
Nevertheless, when considering only the first order of the expansion, the correct result 
for a quantity such as the critical exponent rj can be missed by as much as 60%. Thus, 
as an exp ansion in powers of the field corresponds to an expansion in the number of ver- 



tices 



351 ] . the latter remark is a strong support to approximations, as that of |34j, which 



simultaneously include an infinite number of vertices. 

The LO of the procedure was used in [36| in order to calculate the 2-point function 
of the N = 1 case. Following a simple (and yet accurate) strategy, it was shown that, 
within an analytical and numerical effort similar to that of the DE, one gets a self-energy 
with the correct shape at all momentum regimes: One gets the logarithmic UV behavior; 
its pre-coefficient, which is in fact a 2- loop quantity, follows with only 8% error. In the 
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IR, critical exponents are obtained with a quality similar to that of DE at NLO. As for 
the intermediate crossover regime, a quantity sensitive to this range of momenta was 
calculated to get a result close to the error bars of both lattice and resummed 7-loop 
calculations. 

The purpose of this paper is to apply the method presented in [34|, at its LO, to the 
O(N) model. In (3^] it was shown that the LO of the procedure is already exact in the 
large N limit of the model. Moreover, in this limit, a simple analytical solution of the 
n-point functions at finite momenta was presented. Here, we shall implement the method, 
at any value of N, in order to numerically solve the approximate flow equations of the 
2-point function, at criticality and in d = 3. We shall follow a simple strategy, similar as 
that used in [361 ] . 

The paper is organized as follows. In the next section we shall present the general 
approximation procedure introduced in 34j in the framework of a field theory with TV 
boson fields and, in particular, when the model has O(N) symmetry. In section lITlj we 
shall present two possible strategies to solve the 2-point function equations: the first one 
is simpler, but it looses the above mentioned 2-loop exactness of the effective potential; 
in the second strategy, with a slight increase in the numerical effort, the 2-loop exactness 
is recovered. In section IIVt we present our results, both in the scaling sector and at 
large and intermediate momenta; in particular, we calculate some quantities to gauge 
the quality of the 2-point function thus obtained, and compare our results with those 
following from other means. Finally, in section |V] we study analytically the large N 
behavior of our results, both at leading and next-to- leading order in 1/N, and compare 
them with numerical results of the previous section, as well as with exact results known 
in the literature. 



II. THE APPROXIMATION SCHEME 

In this section we shall briefly present the general formalism of the NPRG and describe 
;he approximation scheme to calculate n-point functions at finite momenta introduced in 



341 ] . We shall make the presentation considering in the first place a generic Euclidean 



field theory with N boson fields (pi, denoted collectively by <p, with action S[(p]. Then, 



4 



we shall specialize to the case where S[<p] has an 0(N) symmetry. 

The NPRG equations relate the bare action to the full effective action. This relation 
is obtained by controlling the magnitude of long wavelength field fluctuations with the 
help of an infrared cut-off, which is implemented I2I, 13 .ll^l I37 ] by adding to the bare 
action S[ip] a regulator of the form 

d d q 



1 f d d q 



(1) 



where (R K )ij(q) denotes a family of "cut-off functions" depending on a parameter k\ above, 
a sum over repeated indices is understood. The role of AS K is to suppress the fluctuations 
with momenta q ^ k, while leaving unaffected the modes with q ^ k. Thus, typically 
(R K )ij(q) ~ n 2 5ij when g < ft, and (R K )ij(q) — > when q 



> 



K. 



One can define an effective average action corresponding to S[(p] + AS K [(p] by r K [</>], 
where <fi is the average field in presence of external sources, <f>i(x) = (ipi(x)). When k — A, 
with A a scale much larger than all other scales in the problem, fluctuations are suppressed 
and Ta[</>] coincides with the classical action. As k decreases, more and more fluctuations 
are taken into account and, as k — > 0, r K=o [0] becomes the usual effective action 



see e.g 



12 



13 



16H . The variation with k of T K [<f>] is governed by the following flow equation 
37J: 

1 f d d q 



rfn 

JJ,l37 



d K TJ 



(27T) 



tr \d K R K (q 2 ) [T^+R K ] 



-1 



(2) 



where denotes the matrix of second derivatives of T K w.r.t. <p (i.e., the matrix of 

5 2 T K /5(j)i5(j)j) and the trace is taken over internal indices. 



components (rl 2 -* 
0: 



For a given value of ft, we define the n-point vertices in a constant external field 



(27r)^ d )(^p,)(rW). 



,12, 



(P 



d d Xl ... / d*x n e i Z3=i p * x i 



(3) 



By differentiating eq. ([2]) with respect to (pi^Xi), ■ ■ ■ , <fii n (x n ) and then letting the field 
be constant, one gets the flow equations for all n-point functions in a constant background 
field cj). These equations can be represented diagrammatically by one loop diagrams with 



\-l 



dt 



1 

'2 





FIG. 1: A diagrammatic representation of the flow equation for the two-point function in an 

represent full propagators and vertices. Crosses represent the 



external field. The lines and dots 
insertion of dtR K 



dressed vertices and propagators (see e.g. 1^]). For instance, the flow of the 2-point 
function in a constant external field reads: 



d d q 



■d t {R K ) in {q) [Gijiq; <p)T { ^ k (p, q, -p - q; <p) 

G k i(q + p; 0)r£i(-p, p + q, -q; 4>)G mn {q\ 4>) 

-^G i:j (q; (j))T { ^ jk (p, -p, q, -q; (j))G kn (q; 0) [ 



(4) 



where G is the matrix of propagators: 



G-'tf-A) = rf(g,-g;0) + R K (q 2 ). 



(5) 



The diagrammatic representation of eq. (jl]) is given in figured! Above, we have introduced 
the dimensionless variable t = ln(ft/A). From now on, as we already did in eq. (jlj), the k 
dependence of the n-point functions shall not be made explicit, unless necessary to avoid 
confusions. 

Flow equations for the n-point functions do not close: for example, in order to solve 
eq. (jJJ) one needs the 3- and the 4-point functions, and respectively. However, in 



ref. 



34| an approximation scheme was introduced in order to solve flow equations yielding 



n-point functions at finite momenta. In doing so, it is possible to exploit two properties 
of these flow equations: i) due to the factor d t R K (q) in the loop integral, the integration 
is dominated by momenta q ^ k; ii) as they are regulated in the IR, n-point vertices 
are smooth functions of momenta. These two properties allow one to make an expansion 
in powers of q 2 /n 2 , independently of the value of the external momenta p. As a typical 
n-point function entering the flow has the form T ( k\pi,p 2 , ...,p n _i + q,p n ~ <?; 4>)i where q 



is the loop momentum, then the leading order (LO) of the approximation scheme consists 
in neglecting the g-dependence of such vertex functions: 

r (n) (Pl,P2, ...,p n -l + q,Pn~ q\ <P) ~ r (n) (p!,p 2 , Pn-1, Vn\ <P) ■ (6) 

Note that this approximation is a priori well justified. Indeed, when all the external 
momenta pi are zero, this kind of approximation is at the basis of DE which, as discussed 
above, turns out to be a good approximation. When the external momenta Pi start to 
grow, the approximation in eq. (El) becomes better and better, and it is trivial when all 
momenta are much larger than k. With this approximation, eq. (HI), for instance, becomes: 

<9jf fe W) = / ^A( R ^(q){G l3 (q-A)r%( P ,0,-p-A) 
xG M (q + p; 4>)T§ ] m (-p, p, 0; <p)G mn {q; 0) 



GiM 4>W%( P , -p, o, 0; 4>)G kn {q\ 4>) \ (7) 



Notice that it is not convenient to also assume q = in the propagators; if this were 
done, the exactness of the LO of the approximation scheme both at one-loop or large N 



limit would be lost (see ref. 



Now, one can exploit the fact that 

<n+l) 



r"' + JI (r> n r> D- th\ - ^i^.-.'"^ 1 ^ 2 ' - Pn ' ] ' , s , 



,(2) 



in order to transform eq. (|7D into a closed equation (recall that G K and T K are related by 
eq. ©): 

*r2(p;*) = f A — ^ L ^^(p,-p; 



^A^M)\G ijM )^^ j 

x G kl (q + P, 0) ^~ P]<t>) G mn {q, 0) 

Eq. (J9j) is valid for an arbitrary theory with iV bosonic fields. In the general case, it 
corresponds to a system of iV(iV + l)/2 equations. From now on, we shall specialize in 
the particular case where the bare action S[ip] has 0(N) symmetry. If one chooses S[ip] 
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to be renormalizable, it is given by 

/I I T %li 1 

\ 2 ( d n*Pi( x ) d ^i( x )) + 2 fi( x )fi( x ) + 4j {<Pi(x)<Pi(?)f \ • ( 10 ) 

In order to preserve the O(N) symmetry all along the flow, it is mandatory to consider 
a regulator respecting the symmetry Doing so, within this approximation scheme, Ward 
identities shall be respected throughout all the flow and, in particular, they shall de valid 
in the (n — > 0) physical limit. The only way to implement this is to consider a diagonal 
regulator. From now on, 

{R K )ij{q) = R K {q)8ij. 

Now, due to the symmetry, the 2-point matrix function can be written in terms of 
only two independent scalar functions; a convenient way to do so is: 

Kb (p> -p; 0; K ) = t a(p; p; «)<U + ^a<PbT B (p; p; «) (11) 

where p(x) = <f) a (x)(j) a (x)/2, a = 1, . . . , N. Following eq. (jSJ), a similar decomposition can 
be done for the propagator matrix. Nevertheless, in this case it proves more convenient 
to use a decomposition in longitudinal and transverse components with respect to the 
external field: 

G ab (p; <p; k) = G T (p; p\ k) ^5 ab - + G L (p; p; «)^^- ( 12 ) 

It is easy to show that 

G T l (p; p; k) =T A (p; p; k) + R K (p), (13) 
G L \p; p; k) =T A (p; p- «) + 2pT B (p; p- «) + R K (p), . (14) 

Using the definition of the functions T A and T B , eq. (TTTT) . as well as that of Gt and 
Gl given above, the flow equation (jUJ) can be decomposed in two equations for Y A and 

d t T A (p; p) = 2pT'li P , p)4% T (P, P) + 2pr|(p; p)4% L (p; p) 

-^((r / A (p;p) + 2 P r>;p))^ I/ ( P ) (15) 

+ ((N- l)T' A (p ]P ) + 2T B ( P ;p))l^ T (p) 



d t r B (p, P ) = (r A (p-,p) + 2r B (p ] p) + 2 P r B (p ] p)) 2 4% L (p-,p) 

+ (N- l)T 2 B (p; p) J® T (p; p) - T'l(p; p) jg T (p; p) 
- r 2 B (p; p) J^ L (p; p) - i ((N - l)T' B (p; p)I® T (p) 
+ (hV' B (p;p) + 2pV' B (p;p))l^ LL (p)) 

+ Tb( P ;p) J ^ d {d t R K (q^B(q;p){G L (q ] p) + G T (q ] p))G L (q ] p)G T (q ] p)} 

(16) 

Above, and from now on, the prime denotes derivative with respect to p and, extending 
definitions already given in 36], we have introduced the functions 

d d q 



An) 
l d;a(3 



(p; «) 



(27r) d 
d d q 
(2n) d 



d t R K (q)GZ- 1 (q;p)G /3 (q;p) 
d t RMG n -\q; P )Gp{p + q-p) 



(17) 
(18) 



with a and (3 standing for either L or T. In eq. f fT6|) . we made use of the identity 



(GUq; P) - G\{q- p)) l - = 2T B {q- p) (G L (q; p) + G T {q- p))G L (q; p)G T (q; p) 

in order to render the expressions manifestly regular at p = 0. 

Eqs. ( ITS]) and ffTB"]) constitute a set of coupled integro-differential equations, with 
respect to the real variables k, p and the modulus of the momentum p. 

Before turning to the strategy to solve it, we shall first comment on an apparent 



inconsistency of this approximation procedure and the way to avoid it [36j. To do so, 
notice that the n-point functions at zero external momenta can all be considered as 
derivatives of a single function, the effective potential V K (p). That is, for example, 

n( 2 )c n N ^ 2 K(p) nn , 

OCpaOCpb 

which entails 

r A (p = 0;p) = — ^— , r B (p = 0;p)= . (20) 

Now, the effective potential satisfies an exact flow equation which can be deduced 
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from that for the effective action, eq. (jSJ), when restricted to constant fields. It reads 

d t V K ( P ) = \J -0y d d t R K {q){{N-l)G T {q-p)+G L { q -p)) = ±{{N-1)I$>. (p; k) 

(21) 

(In the last term of the equation, we made a slight abuse of language with respect to the 
definition of eq. (1T7I) : the function has a unique index a because it contains a unique 
propagator.) 

According to eq. (|T9|) . the second derivative of eq. (}2T]) with respect to the background 
field gives a flow equation for T^(p = 0; p). Now, this equation does not coincide with 
eq. (jSJ) with p = 0: indeed, in contrast to eq. dHJ), the vertices in the equation deduced 
from eq. (12T1) keep all their g-dependence (g being the momentum in the loop integral); in 
other words, it is a more precise equation. There is therefore an apparent inconsistency 
between eqs. (EJ), ( jT9j) and ( 12TT) . However, this can be easily solved. To do so, it is 
convenient to treat separately the zero momentum (p = 0) and the non-zero momentum 
(p 7^ 0) sectors (for a further discussion, see [361] ). 

Let us then define: 

A A (p; p; k) = T A (p; p; k) - p 2 - T A (p = 0; p; k) (22) 
A B (p; p; k) = T B (p; p; k) -T B (p = 0; p; k). (23) 

The flow equations for and easily follow from those for T A and Tb'- 

d t A A {p; p) = d t T A {p; p- k) - d t T A {p = 0; p; k) (24) 

and equivalently for A B (p; p; k). 

The procedure we shall consider in this paper consists then in solving simultaneously 
the 3 flow equations for V K (p) (eq. (1211) ). A A (p; p; k) and A^(p; p; k), and then get 2-point 
functions through 

T A (p; p; «) = p 2 + A A (p; p; «) + (25) 

T B (p; p; k) = A B (p; p; k) + (26) 

The initial conditions for the flow are V\(p) = r^p + (n/6)p 2 (see eq. f lTUj) ) for the 
potential, while those for both A functions are equal to 0, as the classical momentum 
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dependence is explicitly taken out in their definition (see eqs. fl22|) and fl23l) ). Notice 
that, proceeding in this way, not only we maintain the validity of the relationship given 
by eq. (fT9l) . but we also gain more accuracy in the description of the 2-point function: 
indeed, its momentum independent part is now described with a higher precision, while the 
approximation introduced in 34j, which is used in this paper, only affects its j? dependence. 



III. RESOLUTION STRATEGIES 

Although the flow equations for V R (p), A^(p; p; k) and A B (p; p; k) can in principle be 
solved numerically, this nonetheless constitutes a rather cumbersome task. The reason is 
twofold. First, what we called I and J are in fact functionals of the solution r( 2 )(p, p; k); 
this complicates the possible integration strategies. Second, notice that different values 
of p are coupled through the propagators G a (p + q; k) which enter in the calculations of 
the J functions; this demands in principle the simultaneous solution of the equations for 
all p. 

Nevertheless, we shall show that within a simple, and yet accurate, further approxi- 
mation, our flow equations become numerically simpler. In this section we shall in fact 
discuss two possible approximation strategies. 



A. First Level of Approximation: Strategy I 

In order to simplify the above mentioned issues, and to bring down the numerical 
effort necessary to solve our set of NPRG equations, it is possible to perform a further 



approximation (see 36], where the same approximation was used in the N — 1 case and 



an assessment of its accuracy was done). 

Consider first the function I^g(p] k), which does not depend on p. The smoothness 
of the n-point functions and the fact that integrals are dominated by the domain q < k, 
suggest to perform in the propagators of the right-hand-side of eq. (IT?]) an approximation 
similar to that applied to the other n-point functions, i.e., to set q = 0. However, as 
already said in the previous section, in order to maintain both the exact one-loop and large 
N properties of the flow equations, one cannot simply set q = in the whole propagator: 
rather, one needs a momentum dependence recovering that of the free propagators in the 

11 



k —y A limit. Thus, we shall use for the propagators entering the calculation of I^lg{p\ k) 
the following approximate forms 

G T \q; p; k) « Z K q 2 + T A {q = 0; p; k) + R K {q), (27) 
G- L \q- p; k) « Z K g 2 + = 0; p; «) + 2 P r B (g = 0; p; «) + (28) 



where 



dq 1 



(29) 

<J=0,p=po 



As we shall see in the following, the presence of the Z K factor is needed in order to 
preserve scaling properties. As it is well known [14| dT^ 2 '(q; p)/dq 2 \^_ Q depends weakly 
on p. Accordingly, one expects Z K to depend weakly on the value chosen for po. As argued 



in [36( , the choice po = is here the simplest one. With the propagators of eqs. (12711 and 
(J28l) . and the choice of the regulating function [3] 

R K {q) = Z K {K 2 -q 2 )Q{K 2 -q 2 ), (30) 

the function I^a(p; k) can be calculated analytically: 

K d+2-2n f r] K \ 1 

«) = 2 ^^pr- ^ - zf^J (i + ^(p))-(i + ^(p))- (31) 

In this expression, 

r? K = -kO k In Z K (32) 
is the running anomalous dimension and 

rn% M ^-Ai^l = YM (33) 

K Zj k K Zj k 

~2 ( , _ r A (g = 0;p; k) + 2pT B (q = 0; p; K ) V^p) + 2/<(p) ,^ 

m - (p; K) = #z K = ' (34) 

are dimensionless, field-dependent, effective masses. Above, is a number resulting 
from angular integration, K^ 1 = d 2 d ~ 1 Tr d ^ 2 V{d/2) (e.g., K 3 = l/(67r 2 )). Notice that, 

(2) 

for d > 2, I\ a a{K\ p) — > when k — >• 0. 

As for the function J^gip] p; ft), we shall calculate it in a similar way. To do so, let 
us notice that the propagator G a {p + q\ p) in eq. (1T8|) is small as soon as p/tz is large; 
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one can verify that the function J^g(p; p', k) vanishes approximately as K 2 /p 2 for large 
values of p/K. Thus, in the region where J\apip'i P\ K ) has a non-negligible value, one can 
assume p ^ k and then use for G a (p + q; p) an expression similar to that of eqs. (I27j) or 
fl28|) . namely 

G T 1 (p + q; p; k) Z K (p + q) 2 + T A (q = 0; p; k) + i? K (p + g), (35) 
G L \p + g; p; k) « Z K (p + g) 2 + T A (0; p; «) + 2pr B (g = 0; p; k) + R K (p + g), (36) 

One can then calculate the functions </j 3 ^(p; p; ft) analytically (in d — 3 and with the 
regulator of eq. (!30l) ). The resulting expressions are more complicated than those for 
IdapiPi K )i e< l- They are given in appendix [B] Observe that the regulator in eq. fl30|) 
is not analytic at q = k. This generates non analyticities in j]p a p{p\ p; k)] but these occur 
only in the third derivative with respect to p, at p = and at p = Ik (cf., e.g., the odd 
powers of p in eqs. (1B1HB4I) ). and they play no role at the present level of approximation. 

Finally, for the last term in eq. (fT6|) . it is possible to use the same approximation 
procedure described above; doing so, T B (q;p) ~ T B (q = 0; p) = V"(p), and the term is 
then just proportional to the sum of two 1^ functions. 

With the approximation just discussed, both I^lpiP'-, K ) an d jfapiP'^ P'> K ) become ex- 
plicit functions of the potential V K (p) and the field renormalization constant Z K (or, equiv- 
alently, r] K , see eq, (1321) ). From now on, we shall denote the range of momenta p ^ k, 
which is described by quantities as V K (p) and Z K (or r] K ), as the "p = sector" of the 
theory. 

To finish the description of the calculation procedure it is then necessary to make 
explicit how to solve these two flow equations. The simplest way is the usual procedure 
of DE. Specific details are presented in Appendix [A] Here we shall just quote three 
ingredients which are relevant for our present discussion. First notice that, within DE, the 
integral appearing in the potential flow equation is calculated using the approximate 
propagators from eqs. fl27j) and fl28|) ; i.e., the integral is given by eq. fl3T|) . We shall be back 
to this point in the next sub-section. Second, in order to get the proper scaling behavior 
of r^ 2 ) (p; p; k) at small momenta, we need the flow equation for Z K to be consistent with 
the approximate eq. (GO) for the 2-point function. This is achieved by extracting Z K from 
the flow equation of lim p ^ A^(p; p = 0; n)/p 2 , which follows from eq. (Q, and invoking 
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eq. (1291) . Details can be found in appendix [K\ Third, let us notice that, in fact, the 
flow equation for the potential is qualitatively different than those for A^(p; p; k) and 
Ab(p; p; k). Indeed, these depend on a dimensionful quantity p; thus, when k — > 0, the 
corresponding flow stops, giving a finite value for both 2-point functions, which depends 
on p. On the other side, at criticality, the potential only depends on the dimensionful 
physical variable p, whose relevant values shrink to zero when k <C u; the system is then 
characterized by scale invariance. In order to correctly parametrize this property, it is 
convenient to work with dimensionless variables 

p EE K d ^Z^p , v K (p) EE (37) 

Doing so, the dimensionless potential v K (p) approaches a non trivial fixed point form when 

The strategy to solve the eq. ([9]) for the flow of the 2-point function consists then in 
two steps: one first solves the p = sector to get v K (p) and r\ K \ in doing so, the bare mass 
is adjusted in order to reach the IR fixed point. Second, for each value of p, one solves 
the flow equations for the A's, where the kernels I^\p\ k) and Jj(p; p\ k) are explicit 
functions of v K (p), Z K and r] K . The problem of finding the 2-point function of the O(N) 
symmetric scalar field is thus reduced to the solution of a system of partial differential 
equations with parameter p, which can be solved separately for each value of p, and which 
does not involve a numerical effort greater than that required in usual DE calculations. 



B. Improved Approximation: Strategy II 

In section [I] we recalled an interesting property of the approximation scheme intro- 



duced in [34|: when solving the flow equation of the 2-point function at the LO of the 
scheme, the effective potential one gets is exact at 2 loops. Nevertheless, when solving eq. 
(l2Tj) using the propagators described in the previous subsection, this 2-loop exactness is 
lost. We shall present now a simple improvement in the procedure proposed in subsection 
IIII Al in order to recover the 2-loop expression for the potential. This should bring a better 
description of the p = sector of the model, which would be particularly useful for the 
determination of critical exponents. 
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In order to do so, let us briefly remind the origin of the 2-loop exactness. It exploits 
the fact that only one- loop diagrams contribute to the flow (see for example eqs. (jSJ), Q, 
(j2ip . or figured]). Accordingly, the 2-point function gets exact at one loop, provided its 
flow is calculated with quantities which are exact in the classical limit. This is respected 
not only within the approximate eq. ([9]), but also after the extra approximation in the 
propagators introduced in subsection IIII Al As for the potential, eq. (T2"Tj) gives indeed 
an exact expression at 2-loops only if the flow is calculated with quantities exact at 1- 
loop. This is in principle the case for the equations here considered. Nevertheless, the 
extra approximation introduced in subsection IIII Al when applied to the potential flow 



equation, violates this property. Fortunately, this can be swiftly solved, and with a low 
numerical cost. 

This improvement can be achieved by numerically integrating the function 1^ ap- 
pearing in the flow equation for the potential, without using the approximate propagator 
of eqs. (1271) and (|28|) but, instead, the numerical solutions for A^ and A B : 

1 f d d q _ _ . J ... . s 1 



'q* + A A (q;p)+V'(p)+R K (q) 
1 ^ ^ 

+ 



q 2 + A A (q; p) + 2pA B (q; p)+V' + 2pV" + R K (q) . 
One thus needs to simultaneously solve the flow equations for V K , A A and A B . 

Notice that the approximate propagators, and thus the analytic expressions for J, 
from appendix [Bj and for 1^ and I^ 3 \ from eq. fl3T|) . can still be used in the flow 
equations for A^ and A B , without loosing their 1-loop exactness. 

In eq. (|38p . angular integration can be done analytically reducing the problem to a 
numerical integration over one single variable, \q\. Accordingly, the procedure does not 
introduce too much extra complexity in the algorithm. However, an important subtlety 
arises: due to the regulator, the integrand of 1^ takes non-negligible values only in the 
range \q\ < k. Thus, with a fixed grid in q, as k goes to the physical value k = 0, the 
number of points in q for performing a numerical integration in eq. ( |38l) would dwindle 
very rapidly. 

This apparent difficulty is cured by working with fixed values of q/n, i.e., solving the 
A functions flow equations within a grid for fixed values of q = q/n with q < q max ] values 
of q~max ~ 3 — 4 turn out to be large enough. Nevertheless, due to their definition (see 
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eqs. (j22p and fl23|) ). when k — > 0, all functions in the grid would vanish. As usually done 
for the potential (see eq. (15?]) ). this difficulty is simply solved working with dimensionless 
variables: 

Aa(p; p) = ^ , Ab(p; p) = ^ dZ 2 K -i - (39) 

At criticality, these quantities reach finite values in the k — > limit, this limit depending 
on the value of p = p/k. These functions A^ and A# are precisely the quantities entering 
the integrand of the flow equation for the dimensionless potential v K (which follows from 
eqs. (I37j) and (I38D). 

As a final remark concerning the p = sector, notice that we now have many fixed 
point flow equations: those for the dimensionless potential v K and for r] K , and those for 
the dimensionless A functions (in fact, two equations for each value of q on the grid). As 
we are dealing with a Wilson-Fisher fixed point, the flow has only one unstable direction; 
thus once only one bare parameter is fined tuned (here, the bare mass), the complete 
set of equations should reach the fixed point. However, handling the flow equations for 
A^ and A# turns out to be a hard numerical task. For reasons of numerical stability, it 
proves useful to introduce the auxiliary variables 

Ya(p; p) = ~ 2 , Y B (p; p) = . 40 

Flow equations for Ya and Yb are trivially derived. The departing point is the flow 
equations for A^ and A B we used in the previous subsection, i.e., those with the approx- 
imated analytic expressions for the functions J^ 3 \ 1^ and I^ 3 \ Then, making use of eqs. 
(|39|) and (|40|) . a straightforward calculation yields: 

d t Y A (p; p) = V Y A (p; p) + p^{p; ~p) + (d-2 + v )pY A + 2 (l 



x 



Qp V-irj v • ur A • y d + 2 

1 2p w' 2 1 2p w' 2 



p 2 (1 + w + 2pw') 2 1 + w p 2 (1 + w) 2 1 + w + 2pw' 

_ l {Y^p-p) + 2pY'j{p-p)) _ \ {{N-l)Y' A {p ] p) + 2Y B {p ] p)) 
2 (l + w + 2pw') 2 2 (1 + w) 2 

w' 2 

+ 2pJ LT (P; P) (n 2 (P! P)V 2 + 2Y A (p; p)w> + — ) 

w' 2 

+ 2pJ TL (p; p) (Y 2 (p; ~p)p 2 + 2Y B (p; p)w' + — ) (41) 
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Mb,- ~ n 



d t Y B (p; p) = (d-2 + 2 V )Y B (p; p) + p-^(p; p) + (d - 2 + 7/)p^(p; p) 

(N - l)w' 2 (9w 2 + Ylpw'w" + 4pV' 2 ) 



+ 2(1 * 



d + 2 



p 2 (l + w) 3 p 2 (l + w + 2pw') 3 

'2 1 1 „„/2 



11 7i/ 2 1 1 W 



p 2 (1 + w + 2pw') 2 1 + w p 2 (1 + w + 2pw') (1 + wf 
l (N-l) _ 1 (5Y B (p; p) + 2pY B \p; p)) 

2 (1 + w) 2 B[P,P) 2 (l + w + 2pw') 2 

1 Y B (p;p)w' 



1 + w + 2pu>' ^ 1 + y 1 + + 2pw' 

w' 2 

+ (N- 1) J r (p; p) (Kj(p; p)p 2 + p>' + — ) 

r 9«/ 2 
+ J L (p; p)[Y l2 {p- p)f + 4Fj(p; p~)p 2 + 6^(p; p>' + — 

+ 4r;(p; p)r B (pp 2 ; p) + 12Y B (p; p>' + Ap(Y' A (p\ p)Y B (p; p)p 2 

3w'w" 

+ Y A (p; p>" + 3^(p; p>' + + 2F B (p; p)F^(p; p)p 2 

+ 2Y B (p; p>" + 2Y B (p; p)w') + 4p 2 (^ 2 (p; p)) 



it;" 2 

+ 2l£(p; p>" + — ) } - ^lt(p; p) {Y A 2 (p; P)p~ 2 + 2F;(p; p) 

7//2 7// 2 

+ |f) ~ Jtl(P! p) (Yg(p; p)p 2 + 2Y B (p; p>' + ^) (42) 

Notice that the process of going to dimensionless variables brought into play deriva- 
tives with respect to p. Above, we introduced the dimensionless expression J Qj a(p; p) of the 
function J, which is given in appendix [Bj The (k = A) initial values for these functions 
are Y A = 1, Y B = 0. 

In summary, according to this second strategy to get the 2-point functions of the 
O(N) model, one also proceeds in two steps. First, one fixes the p = sector, which now 
demands the simultaneous solution of the full flow equation for the potential (eq. (138]) ) 
together with those for Y A and Y B (eqs. (|4"T1) and (142]) ) for a limited number of values of 
Q < Qmax- Proceeding this way one fine-tunes the bare mass in order to get the infrared 
fixed point. As a second step, one solves the flow equations for dimensionful A a and A B 
(the same one as those of the previous subsection) for any desired value of the external 
momenta p. Within this second strategy, once the p = sector is solved, the p ^ sector 
can still be treated separately for each value of p. 
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Using the alternative idea just presented, one hopes to have a better description of 
the p = sector, and thus, to get an improvement in critical exponents and infrared 
properties of the model. Moreover, one can also hope to retrieve better results at least in 
the cross-over region between IR and UV behavior. 



IV. NUMERICAL RESULTS AND DISCUSSION 



We now turn to the numerical solutions obtained within both of the methods we have 
discussed in the previous section. We shall consider here the properties of the system 
near or at criticality, and we shall study the d = 3 case. Our goal here is twofold: first, 
to asses the quality of the approximation scheme at all ranges of momenta and to verify 
what is the effect of the improvement discussed in the last subsection. Second, since we 
are working here only with the LO of the approximation method proposed in [34] (and, 
moreover, using approximate propagators) we aim to compare our results with those 
already existing in the literature for the O(N) model. 

It is possible to distinguish three momentum regions: the IR sector (p <C u), the 
UV sector (p ^> u) and the intermediate cross-over sector. The first region is dominated 
by scaling properties and has been studied using many different methods with high pre- 
cision. In the next subsection we shall present our predictions for this regime, paying 
particular attention to scaling properties and comparing our results with those existent in 
the literature. As for the UV regime, it can also be described with high precision within 
perturbative calculations. Finally, the crossover region is known with much less accuracy; 
it is therefore the main yield of this work. Both the UV and the crossover regimes shall 
be discussed in the second part of this section. 



A. Scaling properties 

When k C p C «, we expect A^(p; k) to behave as 

p 2 + A A (p; k) = Ap 2 - rf (43) 

where if is the anomalous dimension, i.e., the fixed point value of rj K . As the explicit 
p-dependence of the 2-point function is very hard to obtain, in the literature, the usual 
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way to determine rf is to extract it from the k dependence of Z K (see eq. fl32l) ): 

Z k <xk~v\ (44) 
when k<«. Nevertheless, although these two extractions should in principle lead to the 



same result 



29], when approximating the flow equations this property can be lost. This 



is for example the case in DE which, at criticality, describes physical quantities only at 
p = 0; not even the k <C p <C u regime is correctly described. We have checked that, 
for all the solutions we get, we always have the same result within the two extraction 
methods. 

A more stringent test of scaling is given by the study of the dimensionless quantities 
A A (p; p; K)+p 2 A A (p; p; k) + 2pA B (p; p; k) + p 2 

r> 2 7 ' r> 2 7 ^ ' 

which, in the scaling regime (p, k <C u) should be functions of only p/n and p. In 
fact, according to eqs. fl39|) and PU|) these should be the scaling functions Y A (p; p) and 
Y A (p; p) + 2pY B (p; p). We have numerically checked that our results verify this property. 
In figure [2] we show these scaling functions, as a function of p, for various values of p and 
N = 2. Similar results are found for other values of N. 

Notice that, by definition of Z K (see eq. ( 1291) ). the function Y A (p = 0,p = 0) = 1. 
A non-trivial fact, shown by figure HJ is that both the functions Y A (p,p) and Y A (p, p) + 
2pYb{p v 2) are well approximated by unity for all p < 1 and all relevant values of p (as well 
known [l6|, the latter are those of the order of the minimum of the potential which, within 
our normalization, runs from p ~ iV + 2, when k = A, to p ~ N, when k — *■ 0). This 
is exactly what we assumed within the two numerical approximations introduced in the 

R, the approximated 



361 ]. where the same 



last section. In other words, the figure shows that, even in the deep 
propagators of eqs. (J27j) and (|28j) are indeed accurate. In ref. 
approximation is done in order to solve the N — 1 case, its effect on the functions / and J 
is shown to be small. There, only the longitudinal propagator, i.e., that given by eq. (|28|) . 
plays a role (remember that the term including the transverse propagator is proportional 
to iV — 1). Figure [2] shows that the approximation is even better in the transverse case; 
accordingly, when going to large values of N, one expects the approximations introduced 
in the last section to become better. As we shall see, this is confirmed by our numerical 
results. 
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FIG. 2: The fixed point form for the functions Ya(P,p) and Ya{p,p) + 2pYs(p, p), in d = 3, 
iV = 2, at criticality, as a function of p, for various values of p. 

Our estimates for the anomalous dimension for different values of N, using both of 
our methods of approximation, are presented in table [I] We also plot, in figure [31 Nrj as 
a function of 1/JV. When iV is not too large, iV ^ 4, strategy II introduces noticeable 
improvements. For example, for N = 1, rj* changes from rj* ~ 0.052 using the first 
approximation to rj* ~ 0.047 using the improved one. These values are to be compared 
with results obtained by other means. DE gives 77* = at LO [14] and 77* = 0.033 at 
NNLO [24|. As for the NLO, various results exist: using the regulator of the present paper, 
one gets 77* = 0.050 38], rf — 0.054 with a power-law regulator {22]] and 77* = 0.0467 with 



an exponential regulator [39j; moreover, after an optimization procedure, results move to 
77* = 0.0470, for a theta-like regulator, and 77* = 0.0443, for an exponential one 38]. In 
the table, we present the best results with and without optimization. Thus, even after 
the extra approximation introduced in the numerical resolution strategy, the LO of the 



approximation procedure introduced in 



34| yields a value of the anomalous dimension 



comparable to that of the DE at NLO. All these NPRG values can be compared with 
rj* = 0.034(3) from the resummed 7 loop calculation l|, 77* = 0.0364(2) with high-T 
calculation [3| and rj* = 0.0368(2) from MC fl. 



As A^ increases, our two calculation strategies differ less and less, and eventually they 
give the same result. Most remarkably, our yields are very close to those of the DE at 
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0.025 [40] 
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0.0025 
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0.25/iV 




0.270/iV [411 


1 - 1.034/iV 






1 - 1.081/7V [41] 



TABLE I: Critical exponents for the O(N) model. We present our results within the two 
calculation strategies, together with DE at NLO and the best estimates, with their errors, when 
available. When results of similar quality can be found in the literature, they are both quoted 
(for discussion of DE results, see text). 



NLO, at least for small values of N. Note however that DE results can dep end strongly 



22j depart from 



on the choice of the regulator: with a less reliable regulator, results 
best estimates. Both this work and DE results reach the correct large N limit: rj = 0. 
The large N regime is better studied with the help of figure [3j There, we present the 
best estimates, as well as DE and the present results, together with the known analytical 
numbers for both the large N limit and its first order correction J4I]. One can see that 
both our and DE calculations fail to reproduce the analytical result (15% error within 
DE, 8% for this work). Most remarkably, the large N results seem to be better predicted 
with both NPRG calculation than with the 6-loops resummed perturbation calculation 
of ref. [4fj|; indeed, the latter clearly fails to predict the large N behavior. We shall be 
back to the 1/N limit of our calculation in the last section of the paper. Finally, it is 
interesting to notice that all results present a peak in the value of 77 around iV ~ 2 — 3 
(see tabled]). 

Let us now turn to the critical exponent v. Close to criticality, the effective renormal- 
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FIG. 3: Nrj as a function of 1/N for various model calculations, together with the 1/N analytical 
results. 

ized mass at zero external field (which is the same for the transversal and longitudinal 
propagating modes in the symmetric phase), behaves as 

V'(p = 0) 



|(P=0) 



r — r, 



1 2w 



Ik < u) 



(46) 



where r is the bare square mass and r c its value at criticality. We have performed numerical 
calculations close to criticality in order to determine the critical exponent v for various 
values of N. The results can be seen in table HJ together with those coming from DE at 
NLO and the best estimates in the literature. Within both strategies of calculation one 
gets almost the same results, for all values of N. As occuring with the critical exponent 
T), here comparison of our results with those from DE at NLO depends on the regulator 
used. For example, for N = 1, DE at NLO with the regulator of the present paper gives 



0.625 [381 ]. v = 0.6307 with the exponential regulator ^ 



and v = 0.618 with a power 
law regulator 22|. Optimization makes a tiny difference 38]. In the present work, one 
gets 0.647 (strategy I) and 0.645 (strategy II). These numbers need to be compared with 
0.63012(16) |a| and 0.63020(12) |fj, from high-T and MC calculations, respectively. 

As for the large N behavior, all calculations give the exact limit, v = 1. In order to 
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FIG. 4: N(l — v) as a function of 1/N. This quantity is useful to study the large N behavior 
of the critical exponent v 

strengthen the study of this regime, it is convenient to plot the quantity (1 — v)N, as we 

did in figure HI Notice that these quantity is prone to have much larger numerical error 

than v as one moves to larger values of N. One can see that both DE at NLO and the 

present calculation fails to give the analytical large N prediction for this quantity by not 

more than a few percent. On the other side, as was the case for the other critical exponent, 

I I 

6-loops calculation results from ref. [40| suffer from bigger errors when comparing with 
the exact 1/N expansion. We shall be back to these points in the next section. 

Besides these already very competitive results, the present method is fully devoted 
to calculate quantities that belong to a momentum regime where neither direct pertur- 
bation theory nor scaling properties can be used. These results are presented in the next 
subsection. 



B. Ultraviolet and crossover regimes 

In figure[5]the physical self energy at zero external field, A^(p, p = 0, k — > 0) is plotted 
for various values of N. Notice that a simple dimensional analysis in eq. (jUJ) shows that 
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FIG. 5: A^(p; p = 0)N/(u/N) 2 , in d = 3, for various iV at criticality and zero external field, as 
a function of p/u. As soon as iV ^ 10, the large A limit is attained. 

the self-energy can be written as Aa(p, p = 0, k = 0) = u 2 Aa(p/u, p — 0, k — 0), where 
is a dimensionless function of p/u. On the other hand, in the large iV limit the critical 
self-energy is of order 1/N (see e.g. the next section). As this limit is taken so that uN 
is kept constant, the convenient function to plot is NAa{p/{Nu),p — 0, k — 0)/(Nu) 2 , 
which is shown in the figure. One can see that, as soon as reaches ~ 10, corrections to 
large A^ behavior become small. 

In the perturbative regime (p 3> Nu), one expects r^(p) ~ ((N + 
2)/3)6 a b(u 2 /96n 2 )\og(p/u). In [36| it has been shown, for the A^ = 1 case, that the 
analytical solution of eq. (Q has this behavior, even though the coefficient in front of 
the logarithm, u 2 /9n 4 , is 8% larger (this coefficient is a two- loop quantity and the LO of 
the approximation method introduced in 34j does not include all the 2-loop perturbative 
diagrams exactly). The analysis of [36] is generalizable for all N, provided we include 
the multiplicative factor ((A r + 2)/3)5 a &. Our approximate numerical solutions, within 
both methods, reproduce this result. As explained in [341 ]. at the NLO of the present 
approximation scheme, which is beyond the scope of this paper, the contribution of the 
2-loops diagrams is completely included and the correct prefactor ((N + 2)u 2 5ab/288n 2 ) 
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is recovered. 

For purposes of probing the intermediate momentum region between the IR and the 
UV, we have calculated the quantity 

_ i r d 3 P f 1 1 \ 

u uj (2vr)3 Vp 2 + A A (p) pV 

which is very sensitive to the cross-over regime: the integrand in eq. ( HTl) is peaked at 
p ~ (iV«)/10 |28|. In the 0(2) case and for d = 3, this quantity determines the shift 



of the critical temperature of the weakly repulsive Bose gas [25j. It has thus been much 
studied recently using various methods, even for other values of N. In particular, the large 
N limit for this quantity has been calculated analytically and found to be 1.05 x 1CT 3 



471 ] . As a consequence, 8{(p a (j) a ) /u has been used as a benchmark for non-perturbative 
approximations in the O(N) model. 

In this work, we have found the values for 8(<f> a <f) a ) ju for some representative values of 
N, using both of the methods of approximation for the potential discussed in the previous 
section. The quantity we have chosen to calculate, 8(<f> a <f) a )/u, is the one having a finite 
large N limit. The resulting curves are shown in figure El where they are compared with 



values obtained by lattice calculations [42|, |43j, |44| , resummed 7- loops calculations [45|, |46| 
and to results obtained in ref . 30] . Those numbers, with their corresponding errors when 
available, are also presented in table [III It can be seen that with our approximation 
strategy II one gets slightly better results than with strategy I, but only for small enough 
N. For all values of N where lattice and/or 7- loops resummed calculations exist, our 
results are almost within the error bars of those calculations, except for N = 2, where very 
precise lattice results are available. In the large N limit, one can see that our results differ 
from the exact value by less than 4%. Please note that the large N behavior of 8{<p a <j) a )/u 
is given by 1/N corrections to the self energy 4?j], which are not calculated exactly at 
this level of approximation. Notice also that, within the two numerical strategies we use 
here, one gets the same large N limit for 8{4> a (j) a ) ju. Both in the table and the figure we 



also present results from ref. 30] , corresponding to an iterative ad- hoc calculation within 



NPRG. These results correspond to an improved version of the NLO of that method 



which was a precursor of the procedure presented in [34]. It can be seen that, globally, 
these results are not significantly better than those from the present paper, corresponding 
to an approximate version of the method introduced in [34|, at its LO. 
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7-loops [45J 


ref. [30] 


strategy 


I 


strategy II 


1 


4.94 ± 0.40 x 10" 


" 4 [44] 


4.85 ± 0.45 x 10~ 4 


5.03 x 10" 


-4 


5.44 x 10 


-4 


5.37 x 10~ 4 


2 


5.98 ± 0.09 x 10" 


4 [43] 


5.76 ± 0.45 x 10~ 4 


5.89 x 10- 


-4 


6.44 x 10" 


-4 


6.35 x 10" 4 




5.85 ±0.22 x 10" 


- 4 [42j 
















3 






6.48 ± 0.50 x 10~ 4 


6.57 x 10- 


-4 


7.18 x 10" 


-4 


7.12 x 10~ 4 


4 


7.25 ± 0.45 x 10" 


" 4 [44j 


6.98 ± 0.50 


7.12 > 


< io- 


-4 


7.75 x 10" 


-4 


7.69 x 10~ 4 


10 








8.66 > 


< 10- 


-4 


9.41 x 10 


-4 


9.47 x 10~ 4 


100 














1.10 x 10" 


-3 


1.10 x 10" 3 



TABLE II: Summary of available results for the universal quantity 6(6 a (j) a ) ju. The analytically 



known exact large iV limit for this quantity is equal to 1.05 x 10 



4j| 




FIG. 6: 5{(j) a (j) a ) /u as a function of 1/N 

In summary, when trying to calculate quantities in the momenta crossover region be- 
tween the IR and the UV, the present approximation scheme seems to be particularly well 
suited: already with approximated versions of the LO of the method, one gets numbers 
of about the same quality as those obtained with 7-loops resummed perturbation theory 
or lattice calculations. 
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V. LARGE- iV BEHAVIOR 



As shown in 



34l | . the leading order of the approximating method introduced in that 



paper and used in the present one is exact in the large N limit, for any vertex function. 
This has been verified numerically in the previous section, where the correct large N 
limits, 77 = and v = 1, were extracted from the 2-point function. In this section, in 
order to go beyond these results, we shall study the NLO of the 1/N expansion of the 
model, to calculate both rj and v. The aim is to provide a test for the approximation in a 
situation where it is not exact, as well as to test the quality of the numerical calculations. 

Before turning to this task, it is convenient to remind that the large N limit is taken 
so that uN is constant and then Y K ~ 0(N), V K ~ 0(N), T A ~ 0(1), T B ~ 0(1/N) and 
p ~ 0(N) 41] . In 48| a generic form of the effective action in the large N limit was 
found. Accordingly, a general parametrization of the 2-point function in this limit is 

r$(p, -p; k; <P) = (p 2 + V^p))5 ab + <f) a <p b T B (p, -p; p), (48) 

where r_e is the function defined in eq. ( fill . Using this expression it is possible to show 
that, within this limit, the flow equation for the potential at the leading order of DE (the 
Local Potential Approximation, or LPA) is exact 48j . It has also being exploited to solve 



analytically the NPRG equations for any n-point function 
change of variables [49[] (n,p) — > (k,W), where 

dV 

W = — , 

dp ' 



34j . More precisely, doing the 



(49) 



it is possible to find the 
P = f(W)) Es| and T B 3j: 



UW) =f A (w) + 



N 



arge N limit for both the potential (or its inverse function 
d d q 



u 



F B (p, -p;p) =-= 1 + — 



Nu 



{2k)' 
d d q 



1 



1 



q 2 + W + R A (q 2 ) q 2 + W + R K (q 2 ) 
1 1 



(50) 



6 J (2ir) d q 2 + W + R K (q 2 ) (q + p) 2 + W + R K ((q + p) 2 ) 



-1 



(51) 



which in particular enables for an explicit solution for the f K function (in d — 3 and for 
the regulator ( 1301) ): 



3 3 3 N 
P = f K {W) — -W — -r R + -K d Nn - K d 
u u 2 



2 *TW + \ NK * WC W < 52 > 
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where r# = r\ + KdNAu/3 and C(W) is defined by 



C(W)= dq—±—={ VW ^ 2 , \VwJ J 



Observe that W ~ 0(1). 



A. The critical exponent rj 

The procedure to calculate the critical exponent rj is given by eqs. ( 1All) -( lA6l) of ap- 
pendix [A] From the definitions in eqs. ( lAlj) . ( 1A2[) and ( 1A3I) . one can see that xa ~ 0(1) 
and xb ~ 0(1/N). As for the anomalous dimension, let us first observe that at leading 
order, Xa(p) = 1 (see eq. (148]) ). Accordingly, x^(p) ~ 0(l/iV 2 ). Thus, eq. (1A6[) implies 
that f] ~ 0(1/N). In order to find r/, it is then necessary to get xb at its leading order 
(i.e., 0(l/iV)) and xa at its next-to-leading order (i.e., also 1/N). The first function 
follows from eq. (1A2I) : expanding eq. ( loTT) to order p 2 , one gets: 



Nu 2 K 



(53) 



Finally, using eq. (IA3[) and introducing the dimensionless quantities / = f/(K 3 n), w = 
W/k 2 and u = K^u/n, the exact leading order expression for xb simply follows (remember 
that, for N — ■> oo, Z K = 1 and 77 = 0). 

We now turn to the next-to- leading order expression for xa- In order to do so, we 
shall solve its flow equation, i.e., eq. (1A4I) . up to order 1/N. One can then define: 

1 , *A 
Xa = 1+ AT 

Introducing this expression in the flow equation (1A4|) and keeping only leading terms in 
1/N, one gets: 



dtXA -(d- 2)px A + 



'1 + w) 2 



1 1 + w + 2pw' (1 + w) 2 (1 + w) 2 (1 + w + 2pw') 2 (l + w) 2 K ' 
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Notice that, by construction, xa ~ 0(1) (and then x'a ~ 0(1/N)). Accordingly, in 
the right-hand-side of eq. (I54|) only the leading order expressions of both xb an d w are 
needed. They follow from eqs. fl52l) and fl53l) . after going to dimensionless variables. 

Going to variables (k,w), using <9 t | p = d t \ w + d K wd w and eq. (1A1I) . we end up with 
(where we define d w f K (w) = f' K {w)) 

d t XA{w) - 2wd w XA{w) 

8Nf(w) XB/f'H 2N XB (w) 4Nf/(f~'(w)) 2 



7]N + 



1 + w + 2f(w)/f (1 + w) 2 (1 + w) 2 (l + w + 2/(w)//'H) 2 (l + w) 2 

(55) 



Notice that it is convenient to work in the particular case w — 0, i.e., at the minimum of the 
potential, where the equation is much simpler. Going to the IR fixed point, nd K XA{w) = 
(and, equivalently, u — > oo) it then follows that: 

V* = ^ + o(±) (56) 



AN \N 2 , 

One expects the calculation of 77* to be independent of the renormalization point. 
However, once approximations are done there is normally a small dependence on the 
precise renormalization point. In this case, though, observe that all over this calculation 
we need not to consider the specific renormalization condition introduced in eqs. (1291) and 
(|32|) (i.e., p = 0). It then turns out that the obtained value of Nr/*, when — > 00, does 
not depend on the chosen renormalization point. 

The analytical value of rj* obtained above is reproduced by our numerical results (see 
figure [3] and table 1). The exact value for the 1/N correction to rj* in 0(A) models 



is known 50): 8/(37r 2 ) ~ 0.27. It turns out that the error of our (both numerical and 
analytical) calculation turns out to be of about 8%, which is smaller than the error involved 
in both DE and 6 loops resummed perturbative calculations (see figure [3]). 

A surprising fact is that the ratio between our result (1/4) and the correct number 
(8/(37r 2 )) is exactly the same as the ratio between the correct result for the ultraviolet 
coefficient in front of the logarithm, (A + 2)n 2 /(288vr 2 ), and ours, (A + 2)u 2 /(27tt 4 ). It 



thus seems that, at least for 77, the LO of the approximate procedure introduced in 34| 
misses the second order of the 1/N expansion by the same amount that it misses the 
second order of the perturbative expansion. 
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It is important to notice that the calculation presented above is valid for both of the 
methods of approximation described in section 3 of this paper. This is due to the fact 
that in this section we only needed the large N limit of the potential, which is already 
exact in the case of the approximation used in strategy I. 



B. The critical exponent v 

Let us now move on to the analysis of the 1/N correction for the behavior of the 
critical exponent v. We are dealing here with a critical point for a second order phase 
transition, with only one IR unstable direction under the RG flow. The exponent v is 
related to the eigenvalue of the linearized flow in this relevant direction. 

It will be useful to start with a study of the large N limit for this quantity. Due to 
the significant simplification brought by working close to w = in the study of rj just 
presented, a similar strategy will now be adopted, defining 

b 

w = a{p- po) + g(P - Po) 2 + 0((p - po) 3 ) 



with 



w{p ) = a = w'(po) b = w"(p ) 



It is easy to find equations for this set of variables, which in the large N limit decouple 



481 ]. yielding the equations 



d t po = -(d-2)p + N (57) 

d t a = (d- 4) a + 2Na 2 (58) 

d t b = (2d-6)b + 6N(ab- a 3 ) (59) 

It is then possible to find the non-Gaussian fixed point value for these quantities in the 
large N limit: 

A = ^ 2 (60) 

• (M) 



27V 

6Na* 3 3 (4 - d) 3 



2(d - 3) + 6Na* AN 2 Q-d 
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(62) 



The flow equations for po(t) and a(t), eqs (1571) and fl58|) . can be solved analytically, 
yielding 

p Q {t)=pl + Pe-^ (63) 

°W = a * Ae(^ + l (64) 
These expressions, together with straightforward calculations for the behavior of b(t) 
or any other higher order coupling near the fixed point, clearly show that, in the large N 
limit, the only unstable quantity under the RG flow, as t — > — oo (k — > 0), is po(t). This 
allows for the determination of the expected large N critical exponent v LO = l/(d — 2); 
in particular v = 1 for d = 3. 

Moving on to the analysis of the NLO in the 1/N expansion, we will from now on 
specialize to the simpler case of the first resolution strategy of section IIHI Observe that 
results should not differ by much in the case of the improved strategy of section IIIIBl 
as the numerical calculations for both approximations seem to converge rapidly as N 
increases. 

The NLO flow equation for p Q , for the case of strategy I is 



d t po = -(d-2 + V )p + [ 1 



3 + 2p ^ 

(N-l) + 



(l + 2p a) 2 J 



(65) 



d + 2 / 

Notice here that the functions a(t) and b(t) need only be known in their large iV limit. 
We will focus in solutions for this flow equations that are near the fixed point. We can 
then write, for the quantities in the r.h.s. of eq. (|65l) 

Po(t) = pi + ep^-V (66) 

b(t) = a* + ed(t) (67) 

b{t) = b* + eb(t) (68) 

v (t) =r]* + efip e- t(d - 2) (69) 

where for p and fj we are explicitly taking out the known large N t dependence, and in 
the case of f) we are defining it with an ansatz to be justified later. As for the l.h.s. of 
eq. ( |65l) we can in principle also write po(t) = p* Q + epo^\t)e~ t( - d ^ Expanding eq. ( |65l) 
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to order e one can reach the expression (throwing away terms of order (1/N)) 



9tPo\t) = -7]*p - fip* p - -t^Po + 



d+2' 



o n P t(d-2) 

-Mbe t{d - 2) + Pob* - Plb* a -^—) 
a* a* 



4(3 + ± p *b*)(p* ae« d - 2 » + Po a*) 



1 



;i + 2p*a*) 2 



(70) 



1 + 2p* a* 

As already stated, we know that in the large N limit d(t) and b(t) go to zero after t 
gets negative enough. Nonetheless, we cannot in principle assert the same for the behavior 
of rj, and in fact we will show below that, as rj is already a quantity of order 0(1/N), 
the IR unstable direction of the flow has components in the linearized directions of both 
Po and f]. If we consider \t\ large enough for a and b to be negligible, the following flow 
equation is found 



1 



Po 



'V -VPo 



d + 2 



+ 



2b* 12a* + 8p* b* 
a* 1 + 2p* a* 



1 



,1 + 2pS« 



* n *\2 



(71) 



Knowing the large iV behavior for jr Q (t) (see eq. flTDl ). one can assume that, at NLO 
of the 1/N expansion 



#( t )= A e(-«) 



where, in d = 3, it is easy to check that y\ 



i.e., the correction for v at NLO in 

(i) 



1/N [y = l/(d — 2) + v^-'/N). The equation for p can then be written as an equation 
for 



ZD 



(d-2) 2 N 



-7] -r} Po 



fjN 
d + 2 



2b* _ 12a* + 8p* b* 
a* 1 + 2p* a* 



1 + 2p*a 



* n *\2 



(72) 



expression which, when evaluated in d = 3, using eqs. (1571) . (I6T1) and (|62|) and also (|56|) 
yields the relationship 



v 



(i) 



1 



6Nfj 



(73) 



the —1 in the last expression is the bulk of the required result, as f) is expected to bring 
only a small contribution, which will be calculated in what follows. 

Notice before, that the calculation so far presented is independent of the renormal- 
ization scheme, just as for the case of rj of the last subsection. In what follows we will 
calculate the correction fj, which does depend on the point where the value of r] is fixed. 
Nevertheless, as stated in the last section, the numerical results found for the quantity 
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(1 — u)N suffer from a quite large numerical error for large values of N. This makes 
impossible to distinguish the results between the two strategies presented in this paper, 
or between different renormalization schemes. 

The starting point for the calculation of fj is eq. (1551) for Xa(p, t). In it, scheme 
dependence can only come from the running of r\ K and from boundary conditions. When 
one considers the fixed point solution for xa all the scheme dependence should come from 
the boundary conditions, as we have already shown that rj*, the fixed point solution for 
rj K , is scheme independent. Therefore, the fixed point equation is of the form w8 w xa = 
Nr]* + F*(w), where the r.h.s is scheme independent and thus the only consequence of 
the renormalization scheme on the fixed point solution is an additive constant. 

Near the fixed point, in the r.h.s. F(w,t) = F*(w) + ee"*^ -2 ^ F(w), and similarly 
XA(w,t) = x*a( w ) + te~ ttyd ~' 2) Xa- At order e 

[d - 2)x A {w) + 2wd w x A {w) +Nfj = -F(w) (74) 

Defining the variable xa = Xa + Nr)/(d — 2), it is straightforward to find the only 
solution to the eq. (1741 which is a regular function for all values of w, and also turns out 
to be scheme independent. 

Xa[w) = —-\w\ 2 / aw \w I 2 F(w) (75) 
2 Jo 

As for the function F(w) it is easily obtained by an order e expansion of the function 
F(w), taking into account that (/')* = (/*)', and the fact that xb has a contractive IR 
behavior. 

P( W ) = ~p° n I 8 XbM/(PYM _ i6/*xS»M/((/*)'M) 2 

(1 + w) 2 [ 1 + w + 2f*(w)/(f*)'(w) (l + w + 2f*(w)/(f*)'(w)) 2 
4 1 
" ((/*)'H) 2 {l+w + 2f*(w)/{f*)>{w)f 

16 fW 

((f*)'(w)f (l + w + 2f*{w)/(f*)>(w)y 

where the factor p comes from the t dependence of the function f(w,t) = f*(w) + 
epoe^ d_2 ' ) . Notice that this factor justifies our initial ansatz for f), eq ( 1691) . 

After calculating the integral in eq. (I75|) . which can only be done numerically, one 
can impose the chosen renormalization prescription and obtain the numerical value of f) 
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which would complete the calculation of z/ 1 '. If we take the renormalization condition 
Xa(w(p — 0)) =0 used for the numerical results already presented, we find 

V = ^Xa(w*( P = 0)) (77) 

In d — 3, numerical integration yields xa{ w *{p = 0) ~ 0.028238, which allows for 
finding the NLO 1/N expansion correction for v under our approximation scheme: 

-l-ifi (78) 

which is to be compared with the exact NLO 1/N expansion result [41[ v — 1 — 1.081/iV; 
in this case, then, the analytical error induced by the approximation used in this work 
is of the order of 5%, that is, less than that involved in the determination of rj. Our 
numerical results seem to reproduce the analytical value given by eq. (1781) . even though 
the numerical uncertainty for the quantity N(l — v) grows for N — > oo (see fig. H]). 

As in the case of rj, the results of this work seem to better reproduce the large N 
behavior of the model, when compared with DE or 6-loops resummed perturbative cal- 
culations. In fact, as can be seen in fig. HI there appears to be a notorious disagreement 
between the expected large N behavior and the resummed perturbative results of ref. 

m- 

VI. CONCLUSIONS AND PERSPECTIVES 

In this paper we calculate the 2-point functions of the O(N) model, in all ranges of 
momenta, i.e., either the IR, the UV and the crossover intermediate region. We use a 



method proposed in [34| , which allows for an approximate solution for any n-point function 
NPRG flow equation. Although this method can be systematically improved, in the 
present work we consider only the leading order (LO) of this approximation procedure. In 
fact, in order to deal with a simpler numerical problem, on top of the already approximated 
flow equation we have done extra assumptions to simplify propagators. Moreover, we have 
tested two strategies: while the first one is simpler (we call it strategy I), it misses the 
2-loop exactness of the potential; within the second one (strategy II), which implies a 
small extra numerical effort, the 2-loop expression is recovered. We have shown that both 
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strategies can be justified; in addition, they are exact both perturbativelly and in the 
large N limit. 

We have calculated various quantities to gauge the quality of the 2-point function 
obtained within both strategies. In the IR regime, our solutions correctly reproduce the 
scaling properties. We have calculated two critical exponents, 77 and u, for some represen- 
tative values of N. We got numbers in reasonable agreement with the best known results 
available in the literature, obtained either with lattice or 7-loops resummed perturbative 
calculations. Within strategy II, our numbers turn out to be of the same quality of those 
obtained with Derivative Expansion (DE) at NLO. Notice however that DE is only well 
suited to reproduce these scaling quantities. When going to large values of N, our nu- 
merical results seem to reproduce the analytically known large N limits, exactly at the 
LO and approximately at the NLO in the 1/N expansion. This also happens with DE 



40[ (in fact, 



results but not with those from a 6-loops resummed perturbative calculation 
as soon as iV ^ 5 — 10, results from ref. 40] seem to deviate from the correct large N 
behavior) . 

The UV behavior of the O(N) model 2-point function is of course very well known: 
it presents a logarithmic shape, the pre-coefhcent following from a 2-loop calculation. 
One can analytically prove that the self-energy following from both our strategies does 
have this logarithmic shape. Nevertheless, as the LO of the approximate method of 34| 
does not exactly include all 2-loop diagrams contributing to the 2-point function, the 
pre-coefficent is missed by 8%. We have checked that our numerical results, within both 
methods, reproduce these analytical predictions. 

The intermediate crossover region between IR and UV regimes is certainly the most 
challenging one. In order to check the quality of our 2-point function, we have calculated a 
quantity which is particularly sensitive to this momentum regime. Of course, this number 
cannot be obtained with DE and only lattice or resummed perturbative calculations, as 
well as analytical large N results, exist. In this case, our method is largely competitive 
with the best available results, except for the N = 2 case, where very precise lattice results 
are known. We miss the analytical large N limit for this quantity by 4%; notice that the 
calculated quantity is of NLO in a 1/N expansion, which is not completely included in 
our approximate calculation. 
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Let us finally add a remark concerning our results: as expected, strategy II is more 
accurate than strategy I, but the differences get very small when going to large values of 
N. This is related to the fact that in the large N limit the approximate propagators we 
used in strategy I become exact. 

To our knowledge, this is the first time that the 2-point function of the O(N) model 
is calculated with such overall acceptable properties. Within the NPRG, two ad-hoc 
calculations can be quoted here. First, in refs. [2^] an approximate solution of the flow 
equation including only a limited number of vertices could be calculated, but the result is 
unstable when trying to improve the method 31] . On the other side, in refs 



28 



29 



301 ] . an 



ad-hoc method which can be considered as the precursor of the method presented in 34j 
is used to get a 2-point function with all the expected properties; nevertheless, acceptable 
numbers only follow after an improved NLO order variation of the method, and after a 
lengthy analytical and numerical effort. 

We have also studied analytically the 1/N expansion of the critical exponents 77 and 
v. In fact, as the first order in the expansion (77 = and v = 1) is trivially reproduced 
by our calculation we went up to the second order. In the case of 77, we found a result 
8% smaller than the correct one. As for v, we get (1 — u)N = 1.03 + 0(1/N 2 ), while 
the correct result is (1 — u)N = 1.081 + 0(1/N 2 ). Here, we miss the correct result by 
only 5%. Finally, let us notice that our results seem to reproduce the expected analytical 
predictions. This is a strong support for our numerics. 

To conclude, we observe that this approximate method to calculate 2-point functions 
at finite momenta within the NPRG seems to work properly for all values of N. It 
reproduces, already at the LO, all the expected behavior of the self-energy, giving, in 
most of the regions of N and p, results similar to those of the best accepted values in the 
literature. Moreover, in the large N limit, one gets the best numbers, if comparing with 
DE or resummed 6-loop calculations. The result of this paper thus extend and confirm 



those already obtained in [36] for the Ising universality class. 

For the near future, two works are to be developed. As a natural extension, it would 
be interesting to go beyond strategies I and II and try to numerically solve the full 
approximate equation for the 2-point function, eq. (Q. On the other side, all these works 
call for an application of this method to non-perturbative problems, as e.g. QCD, where 
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only approximate procedures including a limited number of vertices have been considered 
up to now. 



APPENDIX A: THE p = SECTOR 

Within strategy I, one first solves the flow equations of the u p = 0" sector, i.e., those 
for the potential V K (p) and the field renormalization constant Z K (see section IHIj) . In 
fact, this is nothing but a variation of the Local Potential Approximation (LPA), which 
includes a (p independent) field renormalization constant [16]. Although this is a very 
well known procedure, its yields depend on the precise definition one uses for Z K . As 
discussed in the main text, in order to describe correctly the scaling regime, in our case 
the flow equation for Z K has to be compatible with the approximate eq. (Q for the 2-point 
function. 

Let us first consider the flow equation for the potential. It follows from eq. (1211) with 
the functions given by the analytic expression from eq. ( I3T1) . In fact, it is numerically 
preferable to solve the flow equation for w K (p) = dpV K (p) which reads: 

Vk \{(N-1)w' k , 3w' K + 2pw» 



d t w K = -(2-r] K )w K + (d-2 +rj K )pw K 
As for Z K , let us first define: 



d + 2J \(l + w K ) 2 (l + w K + 2pw 



Z a ( P ,k) =1+ — 2 



Zb{p, k) 



dp- 
dA B (p; p) 



p=0 



dp 



2 



(Al) 
(A2) 



p=0 



whose flow equations follow from those of and A#. Then, one imposes the definition 
of Z K given by eq. (1291) . In fact, as the functions Z^(p, k) and Zb(p,k), as well as the 
potential, are zero momentum quantities, it is preferable to define dimensionless quantities 
(see eq. ((371) ): 



Z A (p) ^ ( ^__Z B (P) 

k 4 



Xa(p) =—^—, Xb(p) = . rA _ dZ 2 K -i ■ ( A3 ) 



37 



The flow equations for xa and xb are then: 



dtXA = VXA + (d-2 + r])px' A + 2(1- 
2p 1 



V 



(1 + w + 2pw') 2 l + w 
1 1 



(2X>0 + 



d + 2 
■2p 



2 X bw' 



(1 + w) 2 1 + w + 2pw' 



( x ' A + 2p X 'k)- 



1 1 



2 (1 + w + 2pw') 

2p 1 
(1 + w + 2pw') 2 (l + w) 2 



2 (l + w) 2 



((N-1) X ' a + 2 X b) 



2w 



12 



(A4) 



d tXB = (d - 2 + 2r/) XB + (d - 2 + r/)px' B + 21- 



2(N-1) 

Xbw + 



d + 2 



1 f / / / 

(l + w) d (1 + io + 2pw') 6 I 

+ 4 p(xW + 3xW + 2xz?w") + 8p 2 xW} 

1 2*>' 1 2x B w' 



(1 + w + 2pw') 2 (1 + w) (1 + w) 2 (1 + w + 2pw') 



+ 



1 (W-l) ., 1 

2 (l + w) 2 

1 



2 (1 + w + 2pw') 



—-(5x^ + 2^) 



+ 



1 



1 + w + 2pw' l + w/ l + wl + w + 2pw' 

/2 j 



(N-l) 



w 



(l + w) 4 (l + w + 2pw') 



^— (9w' 2 + 12pw'w A 



+ 4pV' 2 ) + 



2w 



/2 



(1 + w) 2 (1 + w + 2pw') 2 
Finally, the value of ^ K follows from eq. ( 1321) : 

(iV^(p = 0) + 2 XB (p = 0)) 



^2 , (A f X^(p=0)+2 XB (p=0)) 



(A5) 



(A6) 



(1 + W (p = 0)) 2 + 

In strategy I, the p = sector of the theory thus follows from the simultaneous solution 
of the 3 flow equations flATj) . flA4|) and (lA5|) . together with eq. flA6l) . Eq. (lATj) is solved 
starting from the initial condition at k = A: 



w K (p, k = A) = m\ + g A p, 



(A7) 



38 



where the dimensionless parameters m A and g\ are related to the parameters r and u of 
the classical action, eq. (fTU]) . by (d = 3) 

m A = jj, g A = ^-g-, (A8) 

As for the initial conditions for the functions they follow from the definitions in eqs. 
flAT]) . flA2|) . (T52]) and (T23D, together with eq. ([TO]): 

Xa(p;k = A) = 1 , xb(p;« = A) = 0. (A9) 

The parameter r is adjusted in order to be at criticality at zero external field, i.e., 
in order to have a vanishing physical mass. This is achieved imposing the dimensionless 
mass m 2 (p = 0) = w K (p = 0) to reach a finite value when kCm. 



APPENDIX B: THE FUNCTION J 

In this section we present the expressions for the function (p;p;n), defined in 
eq. (1181) . Using the approximate propagators from eqs. (1271) . fl28l) . ( 1351) and ( 1361) the 
integral can be done analytically (in d — 3). One gets (see ref. 36| for the A^ = 1 case): 
a) p > 2,m 2 p < 0. 



1 



kZ 2 k (2tt) 2 (1 + ml) 2 



2 + 



+ p — 3m 



1 

+ 2p 



1 

+ 2^ 



-i + I + (p 



—mi 



2 4 V 



kZ2(2vt) 2 (1 +m2) 2 [p 2 V3 15 



+ 




V_ 

p A \,15 105 



m 



/ + ^)+^(l/(p 6 )) 



(Bl) 
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b) p < 2,m 2 < 0. 



i 



+ 



25?7 
"48 



p 6 + 



7r) 2 (l+m2)2 

4 4n ri 9 

- p + -p + 

3 15 y 3 y 



4 4 



3 77 7r]m^\ 37; _ 2 



12 6 / y 120^ 



2p 
1 



2p 



1 - 



P 



-m 



1 ^ , ^ , 
1 ~2 + A [P 



1 



2 4 




+ 



^(27r)2(l + m2)2 ^3(l+m 2 ) 
2 + 77 — 2m/g + r/m 2 3 2(1 + — 5m^ + r\m 2 ^) 



3(1 + m 2 ) 2 ^ 



6(1 +m 2 ) 3 
c) p > 2,m 2 > 0. 



-P ~ 



tt) 2 (1 + m 2 ) 2 



15(1 + m 2 ) 4 



2 + ^ 



p 4 + 0(p 5 ) 



(B2) 



- + p 2 - 3m J 



-1 + I + (P 2 ~ K) [}-\ + \(P 2 - ™})) - V^y) \ log 

_/ r] 77 . _ 2 2 \ / / m« \ / ma 

-2m@p \1 — — + — (p — rrip) J I Arctan I - — - I — Arctan I — - 



(p-l) 2 +m 2 



(p + l) 2 + m 2 



1 



1 



1 



7/ 



kZ 2 (2tt) 2 (1 + m 2 ) 2 [p 2 \3 15 
+ Y^-i (7 - 35m 2 + V (-1 + 7m 2 )) + 0(l/(p 6 )) 



(B3) 
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d) p < 2, ml > 0. 



(27 



+ m 



2 A2 



' + 4 4 P V 2 8 8 4 



25r)p 3 1 / 4 4?y rip 2 p 3 rjp 3 rjp 5 
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/^ 2 (2vr) 2 (l + <) 2 J 3(1 + m 2 ) 
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f 
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-P 
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Finally, going to dimensionless variables one can define: 



(B4) 



kZI 



(B5) 
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